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Abstract. Suppose a smooth dynamical system has an invariant subspace and a 
parameter that leaves the dynamics in the invariant subspace invariant while changing 
the normal dynamics. Then we say the parameter is a normal parameter, and much 
is understood of how attractors can change with normal parameters. Unfortunately, 
normal parameters do not arise very often in practise. 

We consider the behaviour of attractors near invariant subspaces on varying a 
parameter that does not preserve the dynamics in the invariant subspace but is 
otherwise generic, in a smooth dynamical system. We refer to such a parameter as 
"non-normal" . If there is chaos in the invariant subspace that is not structurally 
stable, this has the effect of "blurring out" blowout bifurcations over a range of 
parameter values that we show can have positive measure in parameter space. 

Associated with such blowout bifurcations are bifurcations to attractors 
displaying a new type of intermittency that is phenomenologically similar to on-off 
intermittency, but where the intersection of the attractor by the invariant subspace 
is larger than a minimal attractor. The presencc of distinct rcpclling and attracting 
invariant sets leads us to refer to this as "in-out" intermittency. Such behaviour 
cannot appear in systems where the transverse dynamics is a skew product over the 
system on the invariant subspace. 

We characterise in-out intermittency in terms of its structure in phase space and 
in terms of invariants of thc dynamics obtained from a Markov model of the attractor. 
This model predicts a scaling of the length of laminar phases that is similar to that 
for on-off intermittency but which has somc differenccs. 

Finally, we discuss some other bifurcation effects associated with non-normal 
parameters, in particular a bifurcation to riddlcd basins. 



1. Introduction 

Nonlincar dynamical systems with invariant subspaces forced by, for examplc, 
symmetries or other constraints are of great physical interest. The dynamics on 
such invariant subspaces can show transverse stability or instability depending on 
whether small perturbations away from the subspace decay or grow with time. As 
these subspaces need not be normally hyperbolic, the dynamics near them can be 
very complicated; for example, the phenomena of riddled basins is typical in such 
systems [Q. 

Recently there have been concerted attempts to understand the bifurcations of 
such attractors on varying a parameter in the system. For system parameters that are 
normal |?J i.e. parameters that leave the system on the invariant manifold unchanged, 
there is a good description of the instability that causes a blowout bifurcation 
; this is a linear instability that can be computed by examining the Lyapunov 
exponents (LE) corresponding to perturbations in transverse directions. The problcm 
of characterising the global branching in such bifurcations is still not understood, but 
some progress has been made in 

In reality, most physically relevant parameters are not normal; that is, they vary 
the dynamics within the invariant subspace as well as that outside it. This, couplcd 
with the fact that most physically relevant chaos is not structurally stable (it is fragile 
in the terminology of ^9)), and the fact that most systems do not possess a skew 
product nature, leads to a variety of novel and physically relevant phcnomcna that we 
investigate in this paper. Thcsc includc 

(i) "Blurring" of a blowout bifurcation caused by the breakup of fragile chaos. 

(ii) Windows within parameter space where a generalised type of on-off intermittency 
(which we call "in-out" intermittency) may appear. 
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(iii) New mechanisms for bifurcation to locally riddled basins, including cases where 
the attractor in the invariant subspace is not chaotic. 

All of these phenomena arise on varying parameters that are not normal in systems 
that are fragile and do not have skew product form. 

More precisely, consider a family of smooth mappings 

/„ : M -> M 

parametrised smoothly by v € R, where M is a compact subset of R m with open 
interior, such that /„ leaves a linear subspace N C M invariant. If the map h — ¡ v \n, 
restricted to N , is independent of v we say that v is a normal parameter for / on N, 
otherwise we say v is a non-normal parameter. 

If a (minimal Milnor) attractor A C N becomes transversely unstable on varying 
a parameter v we call this a blowout bifurcation. For normal parameters this can 
happen at isolated values of v; we will show that for non-normal parameters a more 
typical scenario is the blurred blowout bifurcation, where blowouts can accumulate on 
thcmselves and can even occur on a positive measure subset of parameters. In practice, 
a blurred blowout is recognisable as a complicated pattern of inflating and collapsing 
of the basins of attraction of a family of attractors within N. 

Now suppose that A is a chaotic attractor for / at some particular parameter 
value, and A <£_ N . Suppose moreover that Aq = A n N is non-empty and so thcrc 
are trajectories in A that get arbitrarily close to N but also a finite distance from N, 
arbitrarily many times. In the case that Aq is a (minimal) attractor for /|/v, we say 
that A displays on-off intermittency to the attractor Aq in N (2é| . The attractor A is 
stuck on H to A for /|jv- 

We will see that in general, although Aq is an attractor, it need not be a minimal 
attractor for /|jv; in this case we have "transversely attracting" and "transversely 
repelling" invariant subsets within A ; this more general case we refer to as in-out 
intermittency. If we have an in-out intermittent state that is not on-off, then a minimal 
attractor for /|jv is a proper subset of Aq; these can be, for example, stable periodic 
orbits within N. 

This form of intermittency manifests itself as an attractor where trajectories show 
long periods close to N shadowing orbits in Aq, alternating with short bursting 
phases that may or may not be transient. In particular the "growing" and the 
"decaying" phases can happen via different mechanisms within the invariant subspace; 
in particular, only the phases where the trajectory moves away from N remain close 
to an attractor within N. 

This sort of intermittency is a truly global phenomenon in that it cannot appear 
in systems with skew-product structure, even if they show blowout bifurcations etc. 
Since most systems with invariant subspaces do not have skew-product structure, we 
therefore believe that in-out intermittency will be commonly observable. 

We discuss in Section ^ a simple model mapping of the plane with a non-normal 
parameter. This mapping can be analysed fairly comprehensively, both numerically 
and theoretically. It can be shown to demonstrate all of the above (and many more) 
phenomena. In Section || we put forward arguments to show that the observed 
behaviour of the map in Section ^ is in fact typical. Essential to our arguments are 
the use of Lyapunov exponents (LEs), minimal Milnor attractors, fragility of chaotic 
attractors and the lack of a skew product structure. 
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Section [| discusses and characterises in-out intermittency in more detail, and shows 
that the essential features of this behaviour are well-modelled by a Markov map. This 
map has two parameters that can be used to characterise the intermittency. We dcfine 
a ratio i?¿ of time spent in the "in" and "out" phases. In the limit of i?¿ = we 
regain on-off intcrmittency. This quantity is an invariant of the intermittent dynamics 
up to conjugacy. 

Next, Section ^ examines other bifurcation effects that appear in systems on 
varying non-normal parameters; notably we find a new route to creation of a riddled 
basin via a "non-normal" riddling bifurcation that may well be more common than 
that described in p5fl . 

Finally in Section |^ we discuss our results and point out similar behaviour that 
has previously been found in both ordinary and partial differential equation models. 

2. A model planar mapping with a non-normal parameter 

2.1. The model 

The model we consider extends the well-known logistic map to a mapping of the plane 

f(x, y) = (rx(l -x) + sxy 2 , ue bx y + ay 3 ). (1) 

This has five parameters r £ [0, 4] and (s, v, a, b) £ R and it displays a wide variety 
of bifurcation behaviour. We can view this as a map of R 2 to itself that leaves 
N = R x {0} invariant (in fact for arbitrary parameter values most orbits diverge to 
infinity. Nevertheless, for < r < 4 there is a compact subset of initial conditions 
containing N that remains bounded). The map on N is the well-known logistic map 
and it undergoes the familiar routes to chaos via intermittency and period-doubling 
cascades. 

In the case s = 0, the map has extra structure; it is a skew product over the 
dynamics in x, i.e. it can be written as 

f(x,y) = (h(x),g(x,y)) (2) 

for some g and h, where x € N. We will see that the breaking of the skew product 
form (s t¿ 0) is important to see the generic types of dynamics we report here. 

If we fix r and vary s, v, a and 6 we see that the latter four parameters do not 
affect the map restricted to N; i.e. s, v, a and b are normal parameters for the system 
restricted to N. We are especially interested in the case where these parameters are 
fixed and the only non-normal parameter r is varied. In this case the dynamics in N 
will undergo many bifurcations in regions of interest. 

We investigate the relationship between the dynamics and the numerically 
measured LEs; these are defined as usual for (x, y) <E R 2 and (0,0) ^ (u, v) e R 2 
by 

X( x , y )(u,v) = Hm -\\log DfZ t ) (u,v)\\ 

where this limit exists p3| ; in Section || we will return to explain the dynamical 
behaviour in terms of these rates of asymptotic separation of trajectories. 

Roughly speaking, \ x ,y) ( u -¡ v ) nieasures the exponential rate of growth of 
perturbations in the direction (u, v) along the orbit of (x,y); typically there will be 
two distinct Lyapunov exponents. 
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If (x, y) S N, i.e. y = then we can classify the LE into two cases: the tangential 

LE 

= A( X) o)(tt,0) 

and the transverse LE 

\±(x) = \ x ,o)(u,v) 

for some (u, v) such that u / 0. (Note that it may be the case that \( x .o)( u > v ) — \\\( x ) 
for almost all (u, v); however if it is different for one value of (u, v) this is A^). 

2.2. Numerical experiments 

To get a qualitative understanding of the dynamics possible in non-normal, fragile and 
skew product settings, we numerically calculated a number of dynamical indicators 
for the system (|l|) over a range of control parameters in order to investigate the 
corresponding attractors and their relationship to the invariant subspace N. We now 
summarise some of the most interesting types of new behaviour we have observed. 

2.2.1. Non-normal blowout bifurcation Figure ^ shows the numerically computed 
LEs, for the full system as well as the LE corresponding to the map h(x) = /|jv 
restricted to N and the transverse LE, At, as a function of the parameter r. In the 
top panel we have marked by 1 and 2 two different regions where there are "blowout" 
bifurcations. We show in Figure || typical attractors of system ([[]) with parameter 
values in these regions for some initial conditions. 

We have also calculated in each case the average distances and maximum distance 
of typical trajectories from N and LEs after transients have been allowed to die away. 
These are shown in Figures | and |. 

The regime 1 corresponding to r = 3.828 is brought about by a saddle-node going 
to a stable period 3 orbit in N that breaks down to chaos on reducing r via saddle-nodc 
induced intermittency. The regime 2 located near r = 3.886 seems to have attractors 
that barely vary their statistics. 

In both cases the blowout is observable over a small range of r. As shown in 
Figures [5](d) and ^(d), the transverse LE appears to vary non-differentially over a 
large measure set of r. We suppose that these show essentially the same phenomenon, 
with just a different measure set of points where 'blowout' bifurcations occur. 

Finally, Figure |E| shows the boundaries of the regions in the 2-dimensional (r, v) 
parameter space for which the attractor in N is transversely stable (unstable). Note 
how there are two regions in this parameter space; one in which there is an attractor 
contained in the invariant subspace and another where there is not. The boundary 
between the two takes the form of a graph over r. This is because v is a normal 
parameter and so for a given attractor in N this varies the normal LE smoothly 
through zero. The lack of normality of r means that the variation is much more 
complicated in this direction. 

These numerical computations demonstrate a number of important features of the 
blowout bifurcations observed here, notably: 

(a) There are oscillations through zero of the transverse Lyapunov exponent of the 
attractor for the map h. 

(b) The set of parameter values where the tangential LE is positive (corresponding 
to a chaotic attractor for h) are interrupted by a large number of intervals where 
the LE is negative and the dynamics is periodic. 
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Figure 1. Numerically computed LEs of attractors for the system (gj), upon 
changing the non-normal parameter r. Here Ay and Ai (logistic) are respectively the 
transverse and tangential exponents for an attractor in /|jv (i-e. with initial condition 
y = 0), and Ai and A2 (2D) refer to the LEs for an attractor of the full system. The 
regions marked by 1 and 2 in the top panel correspond to different occurrences of 
"blurrcd blowout bifurcations" , where \t passes through 0. The other parameters 
are fixed at (s, v, a, b) = (-0.3, 1.82, -1, -1) 



(a) 
(b) 



Figure 2. Panel (a) represents the phase space picture for the full system (|lj) 
for parameter values r = 3.8285, v = 1.82, s = —0.3, a = —1 and b = —1 (ín 
the region marked 1), and initial conditions x = 1/2, y = 1/2, whereas panel (b) 
represents a similar picture (corresponding to region marked 2) with parameter values 
r = 3.88615, v = 1.82, s = —0.3, a = —1 and b = —1 and initial conditions x = 1/2, 
y = 1/2. 



Figure 3. Details of the LEs (a,c,d) and the average and maximum distances 
b) of a typical trajectory from the invariant manifold N for attractors of the map 
), in the region marked 1. The minimum distance is always very close to zero 
and hence has been excluded. (d) shows a seemingly non-differentiable change in 
transverse Lyapunov exponent as a function of r, organised by a saddle-node to 
intermittency at r = 3.828. Note however the "spikes" in Ay in the chaotic regime, 
r < 3.282, indicating thc locking onto large periodic orbits that are very close in 
invariant measure to the period three orbit existing for r > 3.282. Increasing the 
number of iterates does not seem to affect the structure of this figure. The othcr 
parameters are fixed at (s,v, a, b) = (—0.3, 1.82, —1, —1). 
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Figure 4. Details of the LEs (a,c,d) and the average and maximum distances (b) 
of a typical trajectory from the invariant manifold for attractors of the map (hj), in 
the region marked 2, where the transients have been allowed to decay. The minimum 
distance is always very close to zero and hence has been excluded, implying that 
all attractors are "stuck on" to N if they are not actually within N. (d) shows an 
apparently non-diffcrcntiable change in transverse Lyapunov exponent as a function 
of r, resulting in a blurred blowout near r = 3.8866. Increasing the number of iterates 
does not seem to affect the structure of this figure. The other parameters are fixed 
at (s, v, a, b) = (-0.3, 1.82, -1, -1). 



Figure 5. Successive scan with amplifications through the two dimensional 
parameter space (v,r), showing regions of transverse stability (below the boundary 
curve) and transverse instability (above the boundary curve) of typical attractors in 
N. The other parameters are fixed at (s,a, b) = (—0.3, —1, —1). This boundary has 
locally the form of a graph over r because v is a normal parameter. 



(c) The minimum distance from N is always observed to be exactly (or very closc 
to) zero, implying that all attractors are "stuck on" to N if they are not actually 
within N. 

On the basis of results in the next section, we conjecture that there is in fact a 
positive measure set of r that correspond to blowout bifurcations in some sense. 

2.2.2. Intermittent behaviour Figure ^ shows two examples of intermittent-type 
behaviour for the system (Q) , involvingan attractor that is not in N but which contains 
points arbitraryclose to N . Figure |6|(a) shows an example of the well-known on-off 
intermittency while Figure [|(b) shows an example of a related intermittency 
(in-out intermittency) where the attractor in N is much smaller than the attractor for 
the full system intersected with N. In the latter case, the statistics of the attractor on 
N are markedly different from the statistics of the attractor for the full system near N. 
This is reflected by the "windows" of periodic behaviour in the x dynamics when y is 
small and growing exponentially over several orders of magnitude. Section ^ presents 
a theoretical analysis and characterisation of this type of intermittent behaviour. 

2.2.3. Non-normal riddling bifurcation Suppose there is a minimal attractor A C N 
with a (locally) riddled basin |lj ||] that is fragile. On varying a non-normal parameter 
it is possible for this to collapse onto a stable periodic orbit within N under arbitrarily 
small perturbations. If the periodic orbit is linearly stable we can infer that its basin 
of attraction has open interior and so upon variation of a non-normal parameter we 
can get transitions from open basins to and from riddled basins in a very natural way. 

This transition to riddled basins will in general be "blurred" as with the blowout 
bifurcation. Figure ^ shows an example of basins of attraction of attractors within 
N computed for the map (0) at two different but close values of r. In one case the 
attractor in N is periodic and in the other a chaotic attractor in N with a riddled basin. 
Between these two values there will be a non-normal riddling bifurcation (or more 
precisely, a set of such bifurcations) . We note that despite their apparent similarity, 
the sets in Figs. R(a) andfj](b) are fundamentally different. To see this we show details 
in Figs. |^(a) and^(b), where the hyperbolic periodic attractor in N clearly shows the 
presence of open sets. 
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Figure 6. Panels (a) and (b) show examples of time series generated by iterating 
(jj) from a randomly chosen initial condition, after transients have been allowed to 
die out. In each case, the top plot shows y in a logarithmic scale, and the second and 
the third plots show x and y on linear scales. Panel (a) shows an example of on-off 
intermittency at r = 3.82786 (Xt = 0.0024) and panel (b) an example of what we 
refer to as in-out intermittency at r = 3.88615 (At = 0.023). The other parameters 
are fixed at (s, v, a, b) = (—0.3, 1.82, — 1, —1). 
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Figure 7. The white regions show the local basin of attraction of (a) a chaotic 
attractor in N and (b) a hypcrbolic periodic attractor in N. In case (a) the basin 
is riddled whereas in (b) the basin is open. Even though the two sets look very 
similar; thcy differ fundamcntally on a finer scale. Points in black arc attracted to 
N without leaving y > 0.1, while those in white cross the thrcshold y > 0.1 before 
going to N. The coordinate region dcpicted by the figures is given by x = [0,1], 
y = [0,0.1], and the paramctcrs are fixed at (r, s,v, a, b) = (3.886,-0.3,0.545,-1,1) 
and (r, s, v, a, b) = (3.88605, —0.3, 0.545, —1, 1) in (a) and (b) respectively. 
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(a) 
(b) 

Figure 8. Dotail of Figurc ^(a) and (b). Panel (b) clearly shows the presence of 
open sets (the existence of the open set can also be inferred by noting that there 
is a linearly stable periodic point at the centre of thc bottom of this figure). The 
coordinate regions depicted here is given by x = [0.49, 0.51], y = [0, 0.05]. 

3. Non-normal parameters: theory 

We explain the new observed effects of non-normal blowout, in-out intermittency and 
non-normal riddling in the context of a recent conjecture of Barreto et al. Q on the 
parameter dependence of chaotic attractors. 

3.1. Fragile chaos and the windows conjecture 

Barreto et al. consider a class of chaotic attractors that are fragile, i.e. attractors 
that are not structurally stable under perturbations, but that nonetheless persist for a 
large measure of nearby parameter values. Such attractors are thought to be present 
in a large number of physically important systems. More precisely, Barreto et al. 
suppose that a mapping h has a chaotic attractor A possessing k > positive LEs. 
Suppose that h' is a mapping that is C 1 close to h and such that almost all points in 
a neighbourhood of A are attracted to periodic orbits. Then they say A is dispelled 
for h'. They define the attractor A as being fragile if A is dispelled for C^-arbitrarily 
close h' . 

In particular, their main conjecture (the windows conjecture) for one parameter 
and one positive LE, suggests that there will be a dense set of nearby parameter values 
at which the attractors are periodic. This set of parameter values is called the windows 
set. 

For thc model (j^) that we consider, there are two deep results about the logistic 
map that mean that the chaotic attractors of the map on N are really fragile in the 
above sense. Namely: 

Theorem 1. (see [jl8|, [Í3|) Consider the mapping f a (x) — ax(l — x) from [0, 1] to itself. 
Thcre is a positive measure subset I c of [0, 1] such that a G I c implies that f a has a 
chaotic attractor with absolutely continuous invariant measure and one positive LE. 
The complement of this set contains an open and dense subset I p C [0, 1] such that 
a 6 I p implics that almost all points are attracted to a stable periodic orbit. 

3.2. Attractors and invariant measures 

We briefly sketch some standard definitions and results that we need. Suppose 
/ : R m — > R m . If A C R m is a compact invariant set then we define 

W S (A) = {x G R m : uj(x) C A} 

to be the stable set or basin of attraction of A; similarly we define 

W U (A) = {x e R m : a(x) C A} 
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to be the unstable set of A. 

Suppose that £ is Lebesgue measure on R™\ Following Milnor p^ j. we say A is an 
attractor if £(W) > 0. It is a minimal attractor if there is no compact invariant proper 
subset A' c A with £(W S (A')) = £(W S (A)). 

Suppose that A is an attractor in the sense of Milnor. We say it possesses a natural 
measure if there is an ergodic invariant measure /j, supported on A such that almost 
all points in W S (A) are generic for fj,; i.e. if almost all x E W S (A) satisfy 

1 í 

hm - V 0(/ fc (x)) = Udfx. 

fc=0 J 

for all compactly supported <fi £ C°(R m , R). 

Note that if A supports chaotic dynamics then typically there are many other 
ergodic invariant measures supported on A\ these are the singular measures: 

M-A = \v ergodic invariant measures : supp(z^) C A} 

and include, for example, measures supported on unstable periodic motions containcd 
within A. 

The basic result of Oseledec's theorem |2Í| is that LEs exist and take a finite 
number of constant values on a set of full í/-measure for any ergodic v. If A has a 
natural measure as well as singular measures we say the attractor is chaotic and thc 
LEs to almost all initial points are given by those for the natural measure; however 
the convergence is non-uniform and shows arbitrarily long deviations towards the LEs 
for any singular measure for A. This causes riddling and several other phenomena 
associated with chaotic behaviour in an invariant subspace. 

As shown, for example, in |[l| pl for a large class of attractors one can 
classify an attractor in an invariant subspace according to its transverse LEs according 
to where zero lies relative to this spectrum. Suppose that \i v is a family of natural 
measures of attractors for h v and dcfinc 

\ T (v) = max{A^(^))} 

the most positive transverse LE for the natural measure at v. Then A v = supp(/z¡,) 
is an attractor for /„ if \ v < 0. It is not an attractor if \t(v) > 0. Parameter values 
vq where \t takes both signs for v arbitrarily close to vq we refer to as blowout points 
of A v and these govern the bifurcation from the invariant subspace [ p6| . 

3.3. Blowout with normal parameters 

The transverse LEs can be thought of as the characteristic exponents for the linear 
skew product system 

f v (x,v) = (h(x),M v (x)v) 
where M(x) is a normal derivative at x; i.e. such that 

D ^ f » ( l ) = ( M v (x) ) ( " ) ' 

where * denotes some matrix function of x. If we consider such a system perturbed by 
normal parameters, the transverse LEs are simply those of a perturbed cocycle and so 
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in particular they can vary continuously with any normal parameter (more precisely, 
they will vary continuously on a generic set of perturbed systems [^j ) . If the LEs do 
vary continuously then in particular Xt will vary continuously and even 

av 

generically for A^ = 0; thus we can see that (at least within an open set of smooth 
systems), generically the set of blowout points has codimension one in parameter 
space. On one side of such a normal blowout, A v is an attractor for /„; on the 
other sidc it is not. When A v loses stability it may give rise to a branch of on-off 
intermittent attractors (non-hysteretic or supercritical scenario), or there may be no 
nearby attractors (hysteretic or subcritical scenario) , as discussed by Ott & Sommerer 
(although there are further possibilities, as discussed in S). Note that the hysteretic 
scenario may give rise to transient on-off intermittent behaviour p4| near blowout. 

3.4- Blowout with non-normal parameters and the blowout set 

If we do not have normal parameters and have fragile chaos in the invariant 
subspace, the attractors A v will typically collapse over some dense set of windows 
and consequently X v will typically not vary continuously. 

As Xt{v) does not vary continuously we cannot apply any intermediate value 
theorem and so it is possible to pass from A^ < to Xt > without passing through 
zero. Therefore, we define 

J+ = {v : X T {v) > 0}, r ={v : X T {v) < 0} 

and the blowout set to bc 

/° = di + ndr = 7+nl~. 

For normal parameters, this is equivalent to the set of blowout bifurcations. The 
numerical evidence of Section || suggests that 1° can be a fractal: more precisely we 
can use Theorem [[] to show the following. 

Corollary 1. There is a map of the form (|l|) whose blowout set 1° has positive Lebesgue 
measure in parameter space. 

Proof: Consider any map that has linear form 

fa(x, y) = {ax(l - x), (a - 2ax)y) 

near the invariant subspace y = 0. This is such that the tangential and the transverse 
LEs are the same for all a. On varying the non-normal a, it is possible to see that 
7+ = I c and I~ = I p in the notation of Theorem [l]; the conclusions of this theorem 
imply that 7° = I + and consequently the blowout set has positive Lebesgue measure. 
QED 

In the more general case of A being an attractor that is stuck on to N with 
Aq = A n N fragile then we conjecture that there are typically positive measure 
blowout sets at loss of transverse stability. If Aq is not fragile but instead structurally 
stable (e.g. if it is uniformly hyperbolic) then the blowout set has zero Lebesgue 
measure. However, even in the case of a positive measure blowout set, the sparcity of 
windows may mean that they are very difficult to observe even in a region where they 
exist. 
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4. In-out intermittency 

We now turn to a static (i.e. non-bifurcation) effect that is associated with the blurred 
blowout. In particular, we characterise a generalised form of on-off intermittency that 
we refer to as in-out intermittency. This includes forms of on-off intermittency whcre 
the attracting dynamics within the invariant subspace does not need to be chaotic. 

In the following, we will suppose we have a minimal Milnor attractor A for / such 
that A = A n N is non-empty and A n N c is also non-empty. 

We write A to be the set of invariant subsets of Aq and define 

Aq = {B C Aq : B is a minimal attractor for /|jv} 

and 

A* = Ubea-B. 

This is a closed invariant subset of the invariant set Aq. We can think of Aq as the 
set of points that attract almost all points within A°. If 

A*q = Aq (3) 

then we say the attractor A displays on-off intermittency. In the more general case 
where 

Aq C A q (4) 

we say that A displays in-out intermittency, implying that the attraction and rcpul- 
sion to the invariant subspace are dominated by different dynamics. 

Hypothesis We will assume that if / has a Milnor attractor A and Aq = A n N is 
non-empty for some invariant manifold N then A is a Milnor attractor for /|j\r- 

Note that if A is an asymptotically stable attractor then this will always hold; at 
present we have no measurable examples where this fails to hold but are not aware of 
a sufficiently general setting where this will always hold. Therefore we will leave it as 
a standing hypothesis. Necessarily this means that Aq will be non-empty, but it may 
be a proper subset of Aq. 

4-1- Skew product systems 

The original examples of on-off intermittency and much theory has been developed for 
skew product systems. Interestingly enough, non-trivial in-out intermittcncy cannot 
occur in such a system, as the following result implies. 

Lemma 1. Suppose that / has skew product form (p|) and A is a minimal Milnor 
attractor for /. Then A = A n N is a minimal attractor for g = f\j>¡. 

Proof: Suppose II : R m — > N defined by II(;e, y) = x is the orthogonal projection 
onto N. If / has form (^) then II o / = g o II shows that g is a factor of /. This means 
that Aq is compact and invariant and moreover, W s (Aq) n N = II(W S (A)); since the 
image of a set with positive m dimensional Lebesgue measure under II has positive 
n dimensional measure, this means that Aq is a Milnor attractor for g. Suppose that 
Aq is not minimal for g. Then there will be a positive measure subset of W s (Aq) n N 



14 



that converge to some proper compact subset B <Z A . Thus there will be a positive 
measure subset of W S (A) that converge to a set B' with H(B') = B and B' C A. 
Taking the closure of B' ', this contradicts the assumption that A is minimal. 
QED 

The following theorem is a direct rcsult of Lemma [l|. 

Theorem 2. If A is a minimal Milnor attractor for a skew product system and A 
displays in-out intermittency then it displays on-off intermittency. 

If / is not a skew product then Aq can be a proper subset of A . For example, see 
the illustration in Figure of a heteroclinic network in R 2 where the only minimal 
attractor in N is a fixed point. 




Figure 9. Heteroclinic cycle consisting of two equilibria in N connections between 
them. If this forms part of an attractor then the only minimal attractor within N is 
the fixed point that is stable within N. 

The existence of in-out intermittency forces several consequences in terms of LEs 
(where At(m) = maxAi(p) for an invariant measure /i): 

Lemma 2. Suppose that A is an attractor displaying in-out intermittency and A^ has 
a natural measure ¡i* for /|jv- Then (i) A is not uniquely ergodic, (ii) At(m*) > 
and (iii) there exists a measure fx' with support on Aq \ A*¡ such that At(a^) < 0. 

Proof: (i) Is a trivial consequence of the fact that Aq ^ A . (ii) Suppose that 
At(a¿*) < 0; then A is a minimal attractor for /. (iii) Likewise, if this were not the 
case then A° would be an attractor for /. 
QED 

As far as the dynamics on A and the minimal attractors A go there seem to be 
many possibilities. In the numerical example, we see examples where Aq is a chaotic 
repellor and A is periodic. There is no reason why A should be composed of a single 
minimal attractor, or indeed why A itself should not be chaotic. 

We now address the question of how to model the in-out intermittent states 
themselves; we do this with the aid of a Markov model. 

4-2. A Markov model of in-out intermittency 

We motivate, describe and analyse a Markov chain model of in-out intcrmittency that 
enables us to predict scalings associated with this form of intermittency. In particular 
we consider the laminar phases, corresponding to phases where the trajectory is below 
a given threshold from the invariant submanifold. Note that unlike the case for on-off 
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intermittency, we do not have a corresponding Fokker-Planck equation which models 
this Markov model (cf |§7)). 

Consider a Markov chain described by two semi-infinite chains of states P¿ and 
Qi, i = 0, 1,2, . . .. We call the chain P¿ the "in" chain and Qi the "out" chain. We 
assign transition probabilities as shown in Table [1] and show the chain schematically 
in Figure [l^. This model has three parameters 

a¡y, flfl and e 

which are all positive and such that a\j + a D + e < 1. For simplicity we assume that 
ajj + a D + e = 1 and so the transition from P¿ to P¿ cannot happen. In this case we 
can reparametrise using A € R by defining 

1 -e + A 1 - e — A 
au = , a D = . 

and note that A corresponds to a bias in the drift on the "in" chain. 
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e a D 
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Pi 
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Qi 








Pi 

Qi 


1 - a D 
1 


- e a D 
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Table 1. Probability of transition from the states P¿ and Q¿ for i > 1 in the 
Markov model of in-out intermittency. The P¿ model a biased random walk that can 
leak onto a deterministic motion away from oo. At i = 1 the Q¿ chain feeds into the 
Pi chain. 




Figure 10. Schematic diagram showing transitions in the Markov model for in- 
out intermittency. The site i can be thought of as indexing a distance f l from the 
invariant manifold for some 7 < 1 and so we approach the invariant subspace on 
letting the index i — > 00. 
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Equilibrium distributions We compute equilibrium probability distributions Pi and 
g¿ of the above model by solving the linear recurrence relations: 

p t = aup i+ i + ot D pi-i + (1 - au ~ a D - e)pi, i > 2 

Qi = qi+i + epi, i > 2 

Pi = &uP2 + qi + (1 - e - a D )p!. 

Expressing in terms of the parameters A and e we use an ansatz pi = // to obtain 

9* = T~~ — M ¿ (5) 
i /¿ 

Pí = (6) 
whcre fx satisfies 

(1 - e + A)/i 2 - 2fi + (1 - e - A) = 0. 
This quadratic equation has solutions 

1 ± VA 2 + 2e - e 2 
1-e + A 

but we are only interested in normalised probability measures; this means that the 
only physically relevant solution (with \fi\ < 1) is 



1 - VA 2 + 2e 

A 4 = 



1-e + A 

For normalisation, the sum over all probabilities is unity 

í = ^C(p t+qi )-C e -^— ±, (7) 

giving C = ^ff^- ■ Now the asymptotic ratio between time spent in thc "in" and 
"out" phases, i?¿ , which is a dynamical invariant, can be easily calculated as 



p _ T, l= i Cpt _ i-/i 

E<=i Cqt e 

and the average value of the transverse variable y can be similarly found to be 

< y >= ¿ C(pí + g0 7 - ,; = ^ 1 ^-- (9) 

This quantity i?¿ can be seen to be an invariant of the system up to conjugacy and 
so is a quantity by which the intermittency can be characterised. The case i?¿ = oo 
(equivalently e = 0) corresponds to when the trajectory stays on the "in" phase and 
hence performs on-off intermittency. In the case i?¿ -C 1, the trajectory will spend 
long times moving away from the invariant subspace on the "out" phase interspersed 
with short and necessarily fast-moving contractions towards the invariant subspace. 
We have seen examples of in-out intermittency in the map ([!]) that show a wide range 
of < Ri < oo. 

Note that log7 corresponds to the transverse LE At for the "out" phases. Usually 
for systems with normal parameters, the average laminar phase as a function of (r— r c ), 
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where r c is the value of r for which the blowout bifurcation occurs, scales as < y > _1 . 
If there is a scaling between < y > or i?¿ and the normal parameters v, et, b or s, 
then we can go further and obtain [i. 

Therefore, we can use this model for any in-out intermittent attractor by 
identifying 

or equivalently, by measuring 



At, A, and e, 



A T , R lo , M and < y > 
and computing 7, A, /i and e via 



/' 

<y> 



1 - VA 2 + 2e - e 2 
1-e + A 

e 



l-/i 
(1-M)7 



7-/x 
At = log^. 

We can also calculate the average length of "out" phases, L a as 



1 — u — e 

i=i r 



(10) 



Note that it would be equally possible to have more than two chains corresponding 
to, for example, more than one type of "out" dynamics. Equally the out dynamics 
could be chaotic and give rise to a biased but two-directional random walk on the 
out-chain. There are many possibilities and we do not examine all of these in detail 
on this occasion, however we have observed such behaviours in the model mapping 
(0). 



4-3. Scalings of laminar phases 

We now compute long time asymptotic properties of the scaling of laminar phases for 
in-out intermittency by analysing the Markov model proposed in the previous section. 
We find that for small e this is very similar to that found for on-off intermittency; 
namely for small times it decays according to a power law n~ 3 ^ 2 whereas for larger 
times it decays much faster, namely exponentially. 

To compute the distribution we consider a trajectory starting at some site on 
the "in" chain and consider the distribution of times as to when it returns for the 
first time to that site in the "in" chain. Suppose we want to find the probability of 
the trajectory returning for the first time after 2n steps. This can occur either by 
the trajectory remaining in the "in" chain for that excursion, or it can happen by 
'leaking' onto the 2n — 2¿th site further down the "out" chain after a time 21 and then 
propagating up the "out" chain for the remaining time. 

More precisely, let 

S^n = P{first return is after time 2n + 2} 
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where this is computed assuming that we start at a given site on the "in" chain and 
Pi n (2n + 2) the probability of the first return occur at time 2n + 2 always in the "in" 
chain. Let 

T2n(2l) = P{attains site 21 after a time 2n, has not returned or switched before then}. 
Then we can compute 

n-l 

S 2n = P in (T=2n+2) + e ^ T 2l (2n-2l). 

l=\n/2~\ 

Using a directed version of the probability of the first return and the hitting time 
formulas g| given by 

P ln (T = 2n + 2) = — j- ( 2n ) (a D au) n+1 
n+1 \ n ) 

and 

r 2n (2/) = ~ ( n 2 ^ y aD r+ i ( au r- 1 

respectively, we obtain, for a¡j 0, 



' ' ' 2 " ^ (a Dau ) n+1 +e J2 V f « ) <P">rW"- ' 



n + 1 V n J — ^ l 

l=\n/2] 



. + uí)(^^fe)" i I ! r(:jw 

It now remains to find the asymptotic behaviour of this expression. Using Stirling's 
formula n\ ~ n" + 2 e^^V^Tr we have 

2?¡ 

n 

which upon taking n — > oo gives 

_3 (4aí/a_D) T 



2/ 



)n i n— 1 _ , 

n n V27rn . 7^ n1 / 
l—\n/2\ 

= h + h- 
We can approximate h by 

/ 2 ^ e f^£V l 1 p ^—^(2xf x +i(2 X -n)- 2x + n -i(au) 2x dx 

\o¿U J n n V2nn J x =\n/2] % 

which by changing the integration variable to x = n(l + y)/2 becomes 

h ~ e f ) ~\ (l + yY^y-y^uf+ydy 

\auj V2Wo y/(l + y)y 

^V^L f\(y)e^dy 
au J V27T J 
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where p(y) = (1 + y) log((l + y)a v ) -ylogy and q(y) = 1 ~ v Note that p(y) has 

\fy\i+y¡ 

a single maximum at y m G [0, 1] given by solving 

p' = log(a¡7(l + y)) - logy = 
and p"(y) = -=\— — |. We need at this point to distinguish between two cases for a¡y. 

Case o/O < ajj < \ In this case, y m = and the maximum y m is in the interior 

of the range of integration. Laplaces' method gives 

f q(y)e n ^ dy ~ q(y m )e n ^ -= 1 

asymptotically asn^oo and so we can estimate 

j e 2au - 1 / a D \" 
2 27r a¡7 — 1 \ 1 — cü¡7 / 

to leading order. Note that this implies that I x <C /2 an d -^2 gives exponential fall off 
of 5*2ti in n for large n. Note however the factor e in /2 means that l\ may dominate 
the statistics up to a moderately large values of n. 

Case o/ 1 > a¡7 > ¿ In this case y m = 1 and p'(l) > 0. This means that the leading 
tcrm in I 2 is given by 

/ 2 ~ en~i (4a D au) n 

as n — > 00. This is the same scaling with n as for /1 and so it will appear to scale 
exactly as for on-off intcrmittency. Thus in this case 

S 2n ~ n-i(4aua D ) n 

in the limit of large n. Rewriting this as 

S 2n ~ n~2 exp(nlog(4a¡7a/j)) 

we note that for any e > 

< 4a¡7a D = (1 - ef - X 2 < 1 
and so S 2n scales as n _3//2 up to the point where exponential decay sets in. 

Other cases In the case a¡7 = \ we note that y m = 1 and p'(l) = 0; we do not 
consider this case in detail. In the case a¡/ = we have a deterministic propagation 
down the "in" chain and it is easy to compute that 

S 2n ~ e(\ - e) n 

i.e. we have exponential decay only, and no n~2 scaling. This case is the limiting case 
where A — ► e — 1. 
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Figure 11. Scaling law for proportion P of laminar phases of length r for in-out 
intermittency at r = 3.8800045 and v = 1.82, a = -1, s = -0.3 and b = -1. This 
case corresponds to a case for which ay < i. Note the n — 3 / 2 fall off changing to 
exponential for larger times. 



4-4- Distribution of laminar phases: summary 

We note that if e > for any value of A the distribution í?2n of laminar phases of 
length 2n always decays exponentially with n. For large n, we consider two cases in 
detail. In the case where ajj > | we have 

S 2n ~ n - 3 / 2 e nlog{4auaD) 

(log(4ar/a£)) = log((l — e) 2 — A 2 ) < with equality if and only if A = e = 0). 

In the case of ajj < i, or equivalently when e > A, there is an asymptotic scaling 
of the form 

y/rr 2tt a\j — 1 

where log(4a¡7Q!£)) < C < 0. This causes the exponential tail to continue to much 
higher n, and indicates the presence of a higher proportion of very long laminar 
phases. The shoulder at high n in Fig. [n] is caused by the contribution of li and 
is a clear indicator of in-out intermittency. Pure on-off intcrmittency, in the noise-free 
case, has a scaling law which is always convex in the log-log plot. It is also possible to 
use these results to understand how the scaling of laminar phases varies on changing 
LEs in the intermittency. 



5. Other bifurcations on varying a non-normal parameter 



The bifurcation to riddling (and associated bifurcation to bubbling) has been discussed 
in detail by |Q |32|, in the case where a normal parameter is varied. 
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In this section we can see how basin riddling can appear on varying a non-normal 
parameter; this occurs via a breakdown of fragile chaos of an attractor within the 
invariant subspace. 

5.1. Bifurcation to riddling 

If we vary a non-normal parameter we expect riddled basins to typically appear and 
disappear in the following way; 
Suppose there is a v such that 

(a) There is a branch of attracting periodic orbits (or more generally, uniquely ergodic 
attractors B v C N for v > v ). 

(b) There is a large measure subset of v < v where there exist a branch of attractors 
C v C N such that B v are in a neighbourhood of the C v and the basins of the C v 
are locally riddled. 

Then we say there is a non-normal bifurcation to riddling at vq. Note that (b) holds if 
there is a transversely unstable measure supported on C v , subject to certain regularity 
conditions. Some of these transversely unstable measures will continue through the 
bifurcation to give rise to complicated repelling sets for v > vq. 

We expect this new bifurcation is in fact a more typical way to create riddled basins 
on varying a non-normal parameter than the normal parameter scenario discussed in 
p5¡ . An example of the basins of attraction before and after such a bifurcation is 
shown in Figure [?]. 

5.2. Other transitions. 

There are many other bifurcations that have been studied in the case of normal 
parameters and which will carry across to non-normal parameter systems, modulo 
the fact that if the chaotic invariant set in N is fragile then these transitions will 
happen on sets that may be considerably large than codimension one. 

For example, the transition to bubbling (3^, [33|, an attractor becoming unstuck 
|| and a transition to cycling chaos || may all appear in a non-normal setting. 

Moreover, the case e — ► in the Markov model for in-out intermittency will 
correspond to a transition from on-off to in-out intermittency. 

6. Discussion 

The motivation for this paper came from previous work of thc authors [|l2) that 
investigated a truncated PDE model of an axisymmetric mean field dynamo model. 
Bifurcation and dynamical effects came to light that did not fit into the usual setting 
of a blowout bifurcation; this led us to consider these effects in the simple mapping 

In this paper we have addressed two main questions; firstly, what happens to a 
blowout bifurcation in non-skew product systems on varying thc dynamics by a non- 
normal parameter within the invariant subspace and secondly, how to characterise 
the dynamics of what we call "in-out" intermittency which we suggest should be a 
commonly observable type of dynamics in realistic models with invariant subspaces. 
We have also raised a number of other questions relating to riddling and other 
bifurcations and discussed these in passing. 
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We have tried to give a general characterisation of in-out intermittency, have 
contrasted it to on-off intermittency and have related its likelihood to the non-skew 
product nature of the system and the non-normal nature of its parameters. We have 
proposed that this type of behaviour can be well modelled by a simple Markov chain 
and have used this to obtain quantitative measures of this state (notably the ratio i?¿ 
of "in"- to "out"- states, which is a dynamical invariant). 

An important signature of in-out intermittency is the presence of intervals with 
exponential growth in the "out" variables with a constant rate, over many orders of 
magnitude. While this occurs, the "in" variables closely shadow the periodic attractor 
in the invariant submanifold, thus giving a clear distinction from on-off intermittency 
where the on and off phases are not so clearly differentiated. 

We shall show in a future publication that in-out intermittency can account for 
the behaviour observed in the study of mean field dynamo models (modelled on partial 
differential equations) and referred to as "icicle" intermittency by Brooke |10(|, who 
uses a skew product model to investigate this (see also Tworkowski et al. |30| and 
Brooke et al. as well as in the truncated versions of related models by the authors 

f¡"2f . Similarly this type of behaviour is related to the "new type of intermittency" 
reported by Hasegawa et al. [[l6| in a ring of phase-locked loops. 

There are still few rigorous results on what has been called "fragile chaos" 
||. Further progress in the rigorous classification of non-normal blowouts is clearly 
hampered by the lack of these results in what is a very difficult area of analysis. 
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